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ABSTRACT 


The  steady  motion  of  an  incompressible  inelastic  non- 
Newtonian  Reiner-Rivlin  fluid  near  a  spinning  cone  has  been 
studied  and  a  similarity  solution  has  been  presented.  It  has 
been  shown  that  the  flow  patterns  can  be  obtained  from 
Srivastava-Jain's  work.  But  the  pressure  distribution  is  not 
the  same  and  we  have  given  the  numerical  values  of  the 
pressure  in  two  tables  and  also  shown  their  variation  in  two 


figures 


THE  BOUNDARY  LAYER  FLOW  OF  A  NON-NEWTONIAN  FLUID 
N^  A  SPINNING  CONE 

Subhendu  K.  Datta 

Introduction .  The  boundary  layer  flows  of  non-Newtonian  fluids  have  drawn 
increasing  attention  in  the  past  few  years  due  to  their  importance  in  many 
technological  fields.  Several  authors  (Srivastava  (1958),  Jain  (1961),  Jones 
(1961)  and  others)  have  solved  the  steady  boundary  layer  flows  near  a 
stagnation  point,  near  a  rotating  disc  or  over  a  plane.  Some  unsteady  boundary 
layer  flows  have  also  been  solved  (Srivastava  (I960),  author  (1961)).  But  to 
my  knowledge  the  three  dimensional  boundary  layer  over  a  spinning  cone  of  a 
non-Newtonian  fluid  has  not  been  studied.  The  problrm  of  three  dimensional 
boundary  layer  near  a  spinning  cone  for  an  ordinary  viscous  liquid  has  been 
studied  by  Wu(1959)«  He  showed  that  with  a  suitable  choice  of  coordinate 
axes  and  Independent  variable  the  flow  functions  can  be  made  to  satisfy  the 
same  set  of  equations  as  were  obtained  previously  by  Cochran-Karman,  but 
the  pressure  does  not  satisfy  the  same  equations.  In  considering  the  same 
problem  in  this  paper  for  a  non-Newtonian  fluid  (incompressible,  inelastic  and 
of  the  Relner-Rivlin  type)  we  have  come  to  the  same  conclusion  that  the  flow 
functions  satisfy  the  same  set  of  equations  which  were  obtained  by  Srivastava 
and  Jain  but  the  pressure  does  not,  when  we  make  a  sulted^le  choice  of  axes. 
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We  have  computed  the  pressure  within  the  boundary  layer  and  shown  the 
variation  in  the  figures  I  and  II. 


Equations  .  The  constitutive  equations  for  an  inelastic  non-Newtonian  Reiner- 
Rivlin  fluid  are 


t;,  =  ' 


l,j  =1,  2,  3 


where  and  e^^  are  the  components  of  stress  and  rate  of  strain  respectively, 

p  is  the  pressure,  and  are  the  coefficients  of  viscosity  and  cross 

c 

viscosity  respectively.  We  know  that 


Figure  1 
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We  choose  cartesian  coordinates  {x,y,2)  with  origin  0  at  the  vertex  of  the 
cone  and  z-axis  along  the  axis  of  the  cone.  Let  a  be  the  semi-opening  angle 
of  the  cone.  Then  if  we  take  coordinates  (ti,  ^,0)  such  that 

X  =  (tj  sinQ;+  £,  cos  a)  cos  0,  y  =  (T^sino  +  4  cosa)  sin0,  z  =ticosa  -  ^sina^  (3) 


then  the  physical  components  of  rate  of  strain  are 

0v  9v,  V  sin«+v„coso; 

e _ a  e - ^  e  -  ^ - 

TjT)  '  8t)  ’  U~  9?.  ’  66  T)Sin«+  4cosa  ’ 


(4) 


,  8v  8v,  ,  8v-  v-sina 

2^84  8t^*’  ^T^0  2  8t|  T|Sino'+  4coso' 


)  , 


.L,!l6  ^0  ^ 

®40  ■  2  ^  84  “  T1  sin  ft  +  4  cos  a 


}, 


where  we  have  taken  the  symmetry  into  consideration. 
Thus  the  equation  of  continuity  is 


8v 


8v 


8ti  84 


V  sin  ft  +  v_  cos  ft 

-a _ L _ 

•q  Sin  ft  +  4  cos  ft 


=  0 


(5) 


The  equation  of  momentum  will  give  after  the  transformation  (3)  , 

2 

8v  8v  V-Sina  „  8t'  8t'  , 

_J1 _ 0 _ j  _  _  .2b  .  Tin  n4 

84  qsinft+  4  cos  ft  '  "  8q  8q  84 


p(v  -r^  +  V. 
T]  8q  4 


qsinft+  4  cos  ft 


,  (6) 
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8v 


P(v 


n  9il 


9t'  dr' 

nL 

C  8£,  Ti sin «+  Ceos  ^  "  8C  8C  8ii 


8v 

Zl. 


V  cos  q; 

O  _ 


cos  Qi(T  -T'  J  +  sin«T' 


HL 


■q  sin  «+  C  eos  a 


(7) 


8v, 


8v, 


9t’  8t' 

\  _  n6  ,  C6 

qsinQ(+  Coosa  8q  8C 


V  sincr+v„cos« 

-a _ _ 


2(sinarT|^Q+cosaT^Q) 
^  qsino+Ccos«  ’ 
(8) 


taking  the  symmetry  into  consideration.  tI^^'s  are  given  by  the  second  equation 
of  (1)  and  equations  (4)  . 

2.  Boundary  Laver  Approximatioa  In  order  to  solve  equations  (5)  -  (8)  we 
assume  that 


v^~0(6/qQ),  Vq'>'0(1),  (9) 

which  are  the  ordinary  boundary  layer  assumptions  (of  interest)  compatible  with 
the  equation  of  continuity  (5)  .  These  assumptions  are  generally  valid  at  a 
large  distance  from  the  vertex  of  the  cone.  Assuming 


0( 


J-  oz 

^0  ^/"sina 


) 


> 


V 

c 


-^OCv)  , 


(10) 


»'  =  l*/p»  =  Pg/p  » 


where  w  is  the  angular  velocity  of  rotation,  we  have  from  equations  (5)  - 
(8)  ,  neglecting  all  terms  smaller  than  0(1)  , 
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9v  dv  V 

'  ^  =  0  , 

3ti  ^  94  Ti 

2  2 

9v  9v  v„  ,  ~  8  V  a  o 

_a  +  — n.+  ^  rX/-JL)2  ^ 

T)  9  T1  4  94  'H  P  ^ 


8,%  !!e  :i6%  ,':a.!lawi.  „!!a>2.,%,2u  ,m 

•%  n  “at* 


''ll  3,,  '^''t  3t  tl  ■■'  J,^2  ■'■■'cUri^St  34>''' 


9v^  9v  V-  9v  9v«  9v  ,  9v-  9v 

X  j-  /_i  .  _Jl  _  _i  .  _J1  _  2—^  •  \  +  -  — - ^  1 

34  ^ai1  94  Tl  94  94  9ti  ^  Tl  94  34  ^  ’ 


1  ^  r  9  >.  ■!  2  .  ^^6 ' 2  ■■  1  .  ti  .2,  .... 

-—  cot<,=  --^  +v^£  — {(-^)  t(-^)  )+cota^(g^  )  ]  .  (14) 


In  deducing  (12)  -  (14)  we  have  used  the  reduced  equation  of  continuity  (11)  . 

So,  now  we  are  to  solve  the  system  of  equations  (11)  -  (14)  subject  to  the 
boundary  conditions 

v^(ti,0)  =  0,  v^(ti,0)  =  0,  Vg(Ti,0)  =  Tiw  sin  a  , 

(15) 

lim  V  (ti,4)  =  0,  lim  v  (ti,  4)  =  0  . 

;-oo  ’I  ;-4.oe  ^ 

In  the  next  section  we  present  a  similar  solution  for  the  system  of  equations 
(11)  -  (14)  subject  to  the  conditions  (15)  . 
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3.  Slmlleurltv  Solution.  We  assume  as  solution  of  equations  (11)  -  (14) 

v^=f(T,)F(C),  v^=g(i,)G(C),  VQ*#(r,)»(0  , 

substitute  these  in  equation  (11)  and  get  as  a  result 


IL&L.  .  _ _ _  A 

GU)  '  ftT,)  +  nf(n)"  2 


(16) 


where  \  is  a  constant. 

To  satisfy  the  third  condition  (15)  we  set  ^t))  s  t]  .  Thus  we  have,  firom 
equations  (12)  -  04)  , 

f(il)  f(Ti)  F^4)  +  g(Ti)  «Iti)  FU)  G(4)  -  “  p  P"U)  + 

+  v^[2f(r,)  f(Ti)  (r(U)^  -  2~{^^^f(0  run  + 

+  ^  (F’(U)^  -  T,^t'(0)^  >]  ,  (17) 

-  nC«(U}^  cot  -  pif  +  ^  n^«’(0)^}+^«Ti))^F'(;))^, 

(18) 

f(ri)  no  t(o  +  3(ti)  0(0  ••(o  +  no  #(0  *  vt,#"(o  + 

+  v^[^(T,^T,)r(0rC0)+^{-2Tif'(T,)r(0n0>  +  2t'(0^n)F'(01  .  (19) 

It  appears  from  the  set  of  equations  (17)  -  (19)  that  the  similarity  solution  is 


possible  if 
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f(n)  *  Ph»C)  *  +  11P2(U  cot  a,  Pj  ~0C1J,  P2  ~0(^)  .  (20) 

Then  from  (16)  we  obtain  g(ti)  =  X  .  Sifstituting  these  in  (17)  -  (19)  we  get, 

2 

equating  coefficients  of  t)  and  t)  , 

+  ^G(4)  F'(0  -  9>\q  =  vF"(y  -  v^[(F'(4))^+2F(4)F'*(4)+3(®'(4))^]  , 


(21) 

~PlU)  *  v^{(F'(4))^  +  (»'(t))^] 

(22) 

^  P^(4)  =  v^(F'(4))^  +  «^(4)  , 

(23) 

2P(&)  «(4)  +  XG(4)  **(4)  =  v«“(4)  +  2v^(F'(4)  »'(4)  ■ 

0 

•  F(4)  9'U)]  .  (24) 

Equations  (21),  (24)  and  (16)  give  us  the  functions  F(4),  G(4)  and  0(4)  , 

while  equations  (22)  and  (23)  give  us  the  pressure  distribution  within  the 

boundary  layer.  Changing  to  dimensionless  variables, 

XG(4)  =  H(4)  , 

H(4)=  Vvw  sin  a  Hj(4^) 
F(4)  »(«sinaFj(4j)  , 

(25) 

«(4)  *  «  Sin  a  »j(4p  , 

P2(^)  =  P"  sin  a  sin  aP2(4j)  . 
1^(0  *  p«^sin^aP^(y,  ^ 
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we  obtain, 

pf (9  +  F-(y  -  (y  =  r-a{l  -  K[(F|(y)^  +  2F^(t^)  FJ-(y  + 

+  3(«;C&j))^]  ,  (26) 

PjC^l)  =  K((FJ(y)^  +  («J(4l))^]  ,  (27) 

P^(y  =(*i(Ci))^  +  K(F'(;j))^  ,  (28) 

2Fj(^)  «j(y  +  Hj(4p  «J(4j)  =  »"i(&j)  +  2KIFJ(4i)  *i(4i)-Fi(y  •"i(Si)l,  (29) 

2Fj(^)  +  H|(4i)  =  0  ,  (30) 

K  s  V  (i>  sin  a/ v  . 
c 

Equations  (26),  (29)  are  the  same  as  obtained  by  Srlvastava  (1958)  for  the 
rotational  motion  of  a  plane  lamina  in  a  non-Newtonian  fluid.  .  The  boundary 
conditions  of  our  present  problem  are: 

Fj(0)  =  0,  »j(0)  =  l,  Hj(0)x0  ,  (31) 

Urn  F  (t.)  X  0,  Urn  *  (t )  x  0  , 

Cj***®  42’*^*  *  * 

and  these  also  correspond  to  the  boundary  conditions  of  his  problem.  The 
expression  for  P2  in  (27)  is  also  same  as  obtained  by  him.  Only  the 
expression  for  in  (28)  differs  from  that  of  his.  This  shows  that  the  nature 
of  the  flow  will  be  the  same  as  in  the  case  of  rotation  of  a  disk.  The  values 

*  Equations  (12),  (13)  . 
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of  Fj,  F|,  and  -H^  have  been  given  in  tables  II-III  by  Jain  (1961)  in  his 
paper  solving  Srivastava's  problem  in  a  different  way.  We  borrow  these  results 
from  his  paper  and  since  the  pressure  distribution  is  different  in  our  case,  so 
we  compute  the  values  of 

Pj/pu^sln^  =  Pj  =  K  , 


and  also  of 


pw 


sina'^v  sin< 


=  p  -  p 
2  20 


0  ^  ^ 


for  K  a  0. 05  and  0. 1,  and  give  these  in  the  following  tables  I  and  n  and  show 
their  variations  in  figures  2  and  3  . 


Table  I 


Pj/p«^  sin^  a 


0.05 

0.1 

0.0 

0.0323 

0.  0646 

0.5 

0. 0160 

0.0320 

1.0 

0. 0081 

0. 0161 

1.5 

0. 0039 

0. 0076 

2.0 

0. 0018 

0.0035 

2.5 

0. 0015 
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Figure  2 


Table  II 


0.05 

0.1 

0.60 

0.4098 

0. 4111 

1.20 

0.5687 

0.  5691 

1.80 

0.6238 

0.6240 

2.40 

0.6418 

0. 6421 

0. 6476 

0. 6479 

3.60 

0.6493 

0.6497 

4.  20 

0. 6499 

0.6502 
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Conclusion.  It  has  been  pointed  out  before  that  the  flow  functions  and 

satisfy  the  same  set  of  equations  as  was  obtained  by  Jain  vdien  he  considered 
the  flow  due  to  a  rotating  plane  lamina  in  a  Reiner-Rivlin  fluid.  So,  the  flow 
pattern  will  be  same  as  in  his  case.  This  means  that  the  values  of  and 

decrease  as  K  increases.  Besides,  the  boundary  layer  thickness  also 
decreases  with  increasing  K  . 

Regarding  the  pressure  distribution  it  may  be  seen  from  table  I  that 
increases  as  K  increases,  but  it  dies  away  vary  rapidly  as  we  go  away  from  the 
wall  of  the  cone.  Also,  from  table  II  we  find  that  (P^  -  P^q)  (^20 
value  of  P^  at  the  wall  of  the  cone)  tends  to  its  asymptotic  value  earlier  as 
K  increases  (which  corroborates  with  the  fact  that  the  boundary  layer  thickness 
decreases).  Table  I  shows  that  P^  dies  away  within  (almost)  half  the  thickness 
of  the  boundary  layer.  Thus  the  effect  of  the  cross  viscisity  is  most  prominent 
only  near  the  wall  and  this  too  becomes  more  magnified  as  we  go  farther  from 
the  vertex  of  the  cone. 

The  numerical  calculations  were  made  in  CDC  1604  of  the  numerical  laboratory 
of  the  Mathematics  Research  Center  by  Jaafar  Al-Abdulla.  My  thanks  are  due 


to  him. 
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